Improving the optomechanical entanglement and cooling by photothermal force 
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Cooling and Entanglement in optomechanical systems coupled through radiation pressure and 
photothermal force is studied. To develop the photothermal model, we derive an expression for 
deformation constant of the force. Exploiting linearized quantum Langevin equations we investigate 
dynamics of such systems. According to our analysis, in addition to separate action of radiation 
pressure and photothermal force, their cross correlation effect plays an important role in dynamics 
of the system. We also achieve an exact relation for the phonon number of the mechanical resonator 
in such systems, and then we derive an analytical expression for it at weak coupling limit. At 
strong coupling regime, we show that utilizing the photothermal pressure makes the ground state 
cooling more approachable. The effect of photothermal force on the optomechanical entanglement 
is investigated in detail. According to our exact numerical and approximate analytical studies, even 
though the photothermal force is naturally a dissipative force, it can improve the optomechanical 
entanglement both quantitatively and qualitatively. 

PACS numbers: 42.50.Wk, 03.67.Mn, 85.85. +j, 42.50.Lc 



I. INTRODUCTION 

Cavity optomecnanics has attracted considerable at- 
tention in the last decade. Such a raising prominence 
is mainly due to its capability in providing new insight 
into the quantum behavior of macroscopic objects and 
exploring the interface between quantum and classical 
mechanics [H-Q. Alongside this fundamental quest, cav- 
ity optomechanics has many other benefits rangingfrom 
very precise measurements on forces and positions |7Hll| 
to quantum information processing [E2h14| . Cooling a 
micromechanical mirror toward its ground state, as was 
achieved at the beginning of this decade [H, El, * s a 
way of reaching the regime where quantum effects are 
important for a macroscopic object. Beyond cooling a 
micomechanical mirror, ways of preparing such a system 
in a nonclassical state have also been suggested [Tj, [l8| ■ 
Furthermore, the entanglement between the optical field 
and the mechanical resonator, something which is of fun- 
damental interest whilst also having important poten- 
tial applications, has also been investigated |19l - l2ll |. The 
fast progress of nanotechnology has enabled the fabrica- 
tion of high quality micromechanical resonators, which 
facilitates cooling of one or more natural modes of such 
resonators to their ground state and entangling the me- 
chanical and optical degrees of freedom of the system. 
Strong coupling between vibrational modes of a mechan- 
ical resonator to optical cavity field, which can lead to sig- 
nificant cooling of the mechanical mode, is demonstrated 
in various types of optomechanical systems [23 - |25| . Such 
a strong coupling is potentially a key feature for creating 
the entanglement. 

From an experimental point of view, a micromechan- 
ical mirror in an optomechanical system can be cooled 
down by two major approaches; self cooling, introduced 
by Braginsky and Vyatchanin [26[, and the cold damp- 
ing. In the latter, an artificial optical feedback is em- 
ployed, which was proposed and investigated theoreti- 
cally by Mancini et al in Ref. [27]. Instead, in a self 



cooling optomechanical device, the off resonant opera- 
tion of the cavity results in a retarded back action on the 
mechanical resonator, and hence, a self modification of 
its dynamics. The optical force which affords this back 
action effect, couples motion of the mechanical oscilla- 
tor to the cavity mode; and such a coupling can become 
strong in appropriate conditions. As mentioned above, 
the strong optomechanical coupling leads to substantial 
cooling of the mechanical oscillator and the optomechani- 
cal entanglement [l9L [28| . The back action is usually gen- 
erated by radiation pressure and photothermal pressure, 
also dubbed in some literatures as bolometric force. De- 
pending on the geometry of the device, photothermal ef- 
fects can become important, e.g. in cantilever suspended 
micromirrors. For the first time it was experimentally 
observed in Ref. [24| that a micromirror mounted on the 
cantilever of an AFM is cooled down due to photother- 
mal effects. An essentially similar device was used in 
later experiments to improve the photothermal cooling 
& 

The photothermal effect was first studied theoretically 
via a classical approach by Metzger et al 30]. In addi- 
tion to the cooling, the self induced oscillation, is also 
observed and studied as the opposite effect [3l|. How- 
ever, cooling a macroscopic or mesoscopic mechanical os- 
cillator toward its ground state brings in quantum effects. 
Therefore, it is necessary to study the problem through a 
quantum mechanical formalism. An interaction Hamilto- 
nian for the radiation pressure effect on a moving mirror 
is provided since 1994 (32l434l |. These works were the ba- 
sis for quantum theory of optomechanical self cooling as 
well as optomechanical instability in the case of radiation 
pressure coupling [35 38]. Nevertheless, a full detailed 
quantum mechanical theory for photothermal coupling is 
still lacking. In fact, in the case of an optomechanical 
system with considerable photothermal effects, because 
of macroscopic degrees of freedom and the dissipative 
nature of the force, the problem is much more sophisti- 
cated. In a semiclassical formalism, the effect was first 
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studied by Pinard and Dantan [39]. They introduced 
a phenomenological Hamiltonian for the photothermal 
force and analyzed the effective mechanical damping rate 
and effective frequency resulting from this Hamiltonian. 
Starting from the same semiclassical model but via a dif- 
ferent approach, the problem is discussed in Ref. [40[ . For 
a setup consisting both radiative and bolometric effects, 
they showed that photothermal force can cool down the 
mechanical oscillator to only five phonons. In an slightly 
modified method it was shown that one can obtain a mi- 
cromechanical mirror nearly cooled to its ground state 

41 1- 

The strong optomechanical coupling regime where the 
optomechanical entanglement can be established, has 
not been studied in the photothermal case, so far. In 
this paper, we start from the phenomenological Hamil- 
tonian of the photothermal effect to derive the quantum 
Langevin equations (QLEs) of the system. Exploiting 
the linearized QLEs, we elaborate a quantum mechan- 
ical analysis for photothermal optomechanics. In fact, 
exact solution of the linearized QLEs provides a formal- 
ism which allows us to describe the system in strong 
coupling regime. To study the optomechanical entangle- 
ment, such a formulation becomes necessary. We shall 
exploit this formalism to study the photothermal cool- 
ing and the photothermal optomechanical entanglement. 
Interestingly, our results state that although photother- 
mal process is naturally a dissipative force, yet exploiting 
the photothermal effect, one can increase value and ro- 
bustness of the optomechanical entanglement. Further- 
more, we present a promising theoretical evidence for the 
photothermal cooling toward the ground state using the 
available devices in the present laboratories. In fact, the 
photothermal force is able to significantly improve the 
cooling, thanks to the strong optomechanical coupling. 

The paper is organized as following: In Sec. II we 
present model of the photothermal force and introduce 
Hamiltonian of the system. In Sec. Ill dynamics of the 
system is studied through standard QLE treatment, then 
we calculate the correlation matrix of the system, and 
from position spectrum of the mechanical resonator ef- 
fective damping and mechanical frequency are derived. 
In Sec. IV photothermal cooling is investigated both an- 
alytically and numerically. The optomechanical entan- 
glement introduced by photothermal pressure is studied 
in Sec. V through numerical and analytical methods. Fi- 
nally, the paper is concluded in Sec. VI. 



II. MODEL 

Essentially an optomechanical system is composed of 
an optical cavity associated with a mechanical degree of 
freedom, which up to now, it has been mounted in sev- 
eral different ways. However, to deal with the problem 
we consider a generic and simple optomechanical model 
consisting a heavy fixed input mirror with a finite trans- 
parency and a light movable mirror constructing a Fabry- 




FIG. 1. (Color online) Scheme of a cantilever suffering pho- 
tothermal effects: the radiation pressure affects on the po- 
sition of the cantilever due to its flexibility, while the pho- 
tothermal force changes the cavity length by the surface ther- 
moelastic waves. 

Perot cavity with a changeable length (see Fig. I). In 
a photothermal system, peculiarly the micromechanical 
mirror is not perfectly reflective, so it absorbs the in- 
tracavity photons. The absorption process leads to ex- 
citation of atomic structure of the mirror, resulting an 
increment in the number of free electrons, which instead 
increases temperature of the mirror. In fact, these pho- 
tons are mostly absorbed at its surface as skin depth of 
the mirror is typically very small compared to the dimen- 
sion of the mirror. The heat diffuses from the surface 
into the mirror through a conduction process. There- 
fore, a temperature gradient develops inside the mirror, 
which is a source for elastic waves (called thermoelastic 
waves) [iH ■ In return, the elastic waves modify length 
of the optical cavity leading to change of the cavity res- 
onance frequency. Therefore, the position of the mirror 
reacts on the cavity field amplitude. The deviation of the 
micromechanical mirror from its equilibrium position is 
a superposition of photothermal modifications, q pt , and 
radiation pressure deformations, q rp [see Fig. I]. 

It is worth mentioning that photothermal effects only 
become important in optomechanical devices in which 
there is significant absorption of the intracavity photons 
by the mechanical oscillator. Strictly speaking, in every 
optomechanical system both photothermal and radiation 
pressure back actions contribute in dynamics of the sys- 
tem. However, the amount of their contribution is dif- 
ferent in various setups. It is practically impossible -or 
at least very difficult- to prepare a setup merely based 
on the radiation pressure or the photothermal pressure. 
Nevertheless, in some schemes the photothermal effects 
are negligible, e.g. an optical cavity with a dielectric 
membrane in the middle [43j , suffers small photothermal 
effects. 



A. Photothermal force 

It is convenient to consider the photothermal problem 
effectively as a feedback force on the movable mirror. In 
fact, it is very difficult to introduce a full quantum me- 
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chanical model for the photothermal process as dealing 
with it needs to consider a large number of degrees of 
freedom; phonons and electrons, and their interaction. 
Therefore, we use a semiclassical model for it. As men- 
tioned above, the effect is a complicated process starting 
from absorption of intracavity photons, leading to mod- 
ification of the cavity field amplitude. Obviously, ampli- 
tude of the cavity field responds to these processes with a 
time lag. Such a delay is dominantly determined by ther- 
mal properties of the absorbing material. These processes 
can be encapsulated in a simple formalism describing the 
action of absorbed photons on the mechanical oscillator 
with an appropriate response function [30| . Therefore, 
the feedback model of the bolometric force is written in 
the following form 



F P Mt)) 



dF dcf (g,s) 
9{t-s) — ds, (1) 



where g(t) = Q(t)(l — exp{— t/rth}) is the response func- 
tion of the system with a thermal characteristic time 
r t h, determined by the heat diffusion procedure. The 
instantaneous deformation force -Fdcf in Eq. ((T|) is phe- 
nomenologically proportional to the absorbed optical en- 
ergy. Therefore, the rate of photon absorption -having 
in mind that the absorption process in the mirror takes 
place with a quantum efficiency equal to j3- is /3/f ut (£) 



with I% ut (t) 



i2 Ut '\t)a% ut (t), where a% ut 



y 2 a - a.j 



is the output field at the micromechanical mirror side 
44]. Taking into account adjustment of the intracavity 
photon energy by position of the mirror, the deformation 
force reads 

F de{ (q,t) = X hLu c {q)m n \t) « X/H«Wc[l - 4^2° Ut W, 

(2) 

where, x is the constant of proportionality and its di- 
mension is inverse of velocity. The deformation constant 
X is discussed in the next subsection where an explicit 
expression for it is derived. 

Integrating by parts, the photothermal force turns 
from Eq. (pj to 

F pt (q(t)) = Fo-J°° d9(t d ~ S) F dcf (q, S )ds (3) 

= ^o + — / &(t-s)e ^ h F dei (q,s)ds, 

r th J-oo 

where Fq is a constant force whose only effect is modifi- 
cation of the mechanical equilibrium position. 



B. Deformation constant 

Let us discuss about the deformation constant intro- 
duced in the photothermal force equation @ . In Ref . [39| 
through a phenomenological assumption the constant 
was taken equal to 2/c, where c is the speed of light. 



Instead, in Ref. [4l|, the authors argued that the con- 
stant considered in Ref. [39| does not provide the actual 
strength of the force. Therefore, they introduced a factor 
which pertains to the quantum efficiency of the photon 
absorption (in our notation it is /3) yet it can become 
much greater than unity. On the other hand, Dc Libcrto 
et al obtained a value for x by fitting their theory to 
the experimental data [Io|. Here, we shall suppose that 
/3 gets values between and 1, and then we show that 
X essentially relates to elasticity of the micromechanical 
mirror. Starting from the thermoelastic wave equation, 
we find an explicit expression for %. In the following 
lines, we temporarily suppress radiation pressure effects 
and focus only on the photothermal force. 

The three dimensional thermoelastic wave equation 
can be split into longitudinal and transversal components 
respect to its direction of propagation (see appendix A) . 
Here we are only concerned about the longitudinal com- 
ponent, since only it is responsible for modifications of 
the cavity length. Therefore, assuming x as symmetry 
axis of the cavity, the thermoelastic wave equation is 
given by 



dx 2 



E Y ath dT 
3g(I - 2o)~dx~ 1 



(4) 



where Ey is Young modulus, a is Poisson ratio, ath is 
linear heat expansion coefficient of the micromechanical 
mirror and g is its density. Note that the source term 
on the right hand side of Eq. (j4]) involves the tempera- 
ture gradient inside the mirror which in turn stems from 
the heat diffusion. For the sake of simplicity, we have 
assumed that there is a temperature gradient only along 
the cavity axis. This is reasonable for a uniform heat flux 
on the surface of the mirror. Since surface of the mirror 
is heated by the photon absorption process, considering 
a uniform cavity field profile makes our supposition reli- 
able. In fact, one can define an effective radius for the 
cavity beam profile (which is typically a Gaussian) which 
its prevailing power appears inside a circle with radius rg . 
Thus, this assumption is valid for a mirror with area small 
compared to cross sectional area of the cavity beam. 

On the other hand, neglecting effect of deformations on 
the temperature of the material, thermal diffusion equa- 
tion is given by 



gC 



dT 
~dt 



A", 



i.h 



d 2 T 
dx 2 



Q 



tin 



(5) 



where C is specific heat capacity of the material, K t h is 
its thermal conductivity and Qth is the heat flux. The 
heat flux is given by 



hh{x,t) =^7 2 n eo|E(t)| 2 (5 s d<5[a;- 



«(«)]■ 



(G) 



where no is the refraction index and <5 s( j = (2//j,oc e iWL) 1 ^ 2 
is the skin depth of the material for frequency of the 
laser wl- Also, a c \ is electrical conductivity of the mir- 
ror, and E(i) is the classical intracavity electric field. In 
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photothermal setups the surface of the cantilever is usu- 
ally coated by a thin layer of gold (see e.g. [Hj]) which 
its penetration depth for optical frequencies is about a 
few nanometers. Hence, it is sensible to assume that the 
cavity photons are absorbed just on the surface of the 
mirror. This assumption is imposed through the Dirac 
delta function in Eq. ([6]). 

From ([5]) and (JSJ) an expression for the temperature 
gradient can be derived, and then replaced in the right 
hand side of the elastic wave equation (j4|). To perform 
Eq. Q we use spatial Fourier transform of the longitu- 
dinal deformation u T , defined as 



u x (x,t)= I dk x e tk * x u x (k x ,t) 



(7) 



It is clear that the cavity length is only modified by am- 
plitude of the elastic waves on the surface of the mirror, 
therefore, one sets x = to arrive at 



dk 



X ^J^X 



3gK th (l - 2a) 



Wl 2 } 



-•;», ~ — : — n eo|E(i)| 2 \ = 0, 

where we have introduced the longitudinal elastic wave 
frequency as uj x = k x c x . Now let us define the position 
of the micromechanical mirror q(t) as superposition of 
elastic wave components on the surface, that is 



?(*) 



u x (k x , t)dk x . 



(9) 



Therefore, including the fluctuation and dissipation pro- 
cesses, from Eq. ([5]) one arrives at the following equation 
of motion 

q + j x q + u 2 x q = — (J"th + -7 r P t), (10) 

TO 

where J^h is a classical thermal Langevin force and 



3K th {l - 2a) 



- 72 n eo|E(tr, (11) 



is a classical interpretation for the photothermal force, 
whose time lag is implied in |E(i)| 2 (see Ref. [45j]). A 
comparison between (fTTj) and (|2|) leads us to the following 
relation for the deformation constant 



3tf t h(l - la) 



(12) 



We consider the photothermal setup in Ref. [29j to es- 
timate the value of x- The cantilever utilized in Ref. [29| 
is composed of a silicon substrate with a few nanometer 
gold coating. This gold layer is too thin to contribute 
significantly in elastic vibrations of the cantilever. Ac- 
tually, the temperature gradient effectively ranges from 
the surface of the mirror through the thermal diffusion 
length Z t h = (ZnKth/ qCuj-l) 1 / 2 . One can easily estimate 
that the thermal diffusion length of gold is a few microm- 



eters, while as mentioned above the gold coating on the 
cantilever is just a few nanometers. Thus, the silicon 
substrate dominantly contributes in the thermal expan- 
sion (and contraction) of the mirror. Inserting thermal 
and elastic parameter values of silicon into Eq. (fT2"|) one 
gets x ~ 10~ 5 , which is almost the value obtained by the 
curve fitting of Ref. (4p| | and considered phenomenologi- 
cally in Ref. [H|. 



C. Hamiltonian of the system 

As discussed above, in an optomechanical system the 
effect of the bolometric force can be realized as a feedback 
process. Therefore, the effective photothermal Hamilto- 
nian can be written from Eq. @ as 



Hpt(t) = -_ q{t) 



nh 



h(t-a)[l-^f-]l^(a)da, 

where w c x/?/rth functions as the feedback loop strength 
and /2 Ut ( s ) as the corresponding output photocurrent. 
In practice, this photocurrent is absorbed by the mir- 
ror and acts on momentum of the mechanical resonator. 
Moreover, the response function of the system from ^ 
is h(t) — O(t) exp{—t/r t h}. To investigate quantum dy- 
namics of the system, we must write the Hamiltonian of 
the whole system. We consider an optomechanical sys- 
tem composed of a mechanical resonator with frequency 
w m which is subject to radiation pressure and photother- 
mal pressure due to an optical cavity mode with fre- 
quency uj c which is driven by an intense laser. It is conve- 
nient to transform canonical operators of the mechanical 
oscillator so that [q,p] — i. Hence, the Hamiltonian reads 

H = fiiu c a) a + —hu m (p 2 + q 2 ) — KG^a) aq 



HXq / h(t - s) 



G Q q(s) 



l2 Ut (s)d S 



(14) 



where the first term describes the cavity mode, with an- 
nihilation operator a ([a,a^] = 1), while the second term 
corresponds to the free mechanical oscillator. Third term 
of (|T4"]) is the radiation pressure interaction, while the 
fourth represents the photothermal Hamiltonian. The 
last one is the effective Hamiltonian of a input laser 
with power P and frequency ojl- In Eq. (fl4"|) . Go = 
{uj c / L^^Jh/rauj^ is the single photon radiation pressure 

coupling and A = (x/3w c /r t h) \/ W/ muj m = x/3LG /T th is 
the photothermal coupling strength, while the mean am- 
plitude of the input field is related to the laser power by 
£ = \J2^\P '/huL, where 71 is the cavity loss rate through 
its input port. We notice that since typical optomechan- 
ical parameters suggest that Goq <C U) c , we can neglect 
the second term of integrand in the photothermal Hamil- 
tonian. 
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III. SYSTEM DYNAMICS 

A. Quantum Langevin equations 

A full description for dynamics of the optomechanical 
system is provided by its QLEs which can be derived 
from Hamiltonian (1141) . Fluctuation and dissipation pro- 
cesses can be taken into account in a fully consistent way 
[46[ . We write these QLEs in rotating frame of the laser 
frequency as the following 

q = uj m p, (15a) 
p = -oj m q - 7,„p + G a^a 

+ ^G / h(t - s)I° ut ( S )ds + {, (15b) 

7th J 

a = -(7 C + iA Q )a + iG qa + £ + \[ 7 2q~\ a 1 " 

+ v /2^ 4 n , (15c) 

where Ao = cj c — wl is detuning of the laser frequency 
from the corresponding cavity mode, while 71 and 72 are 
the cavity decay rates from the heavy input mirror and 
the micromechanical mirror, respectively, constructing 
the entire cavity decay rate 7 C = 71 + 72 ■ Furthermore, 
£, a™, and a™ are the zero mean value Langevin noise 
operators affecting the system: The mechanical mode is 
affected by a Brownian stochastic force £ it) , obeying the 
correlation function 47 1 



mm) = ^ / ^< 

w m J 2n 



1 y 2k B T Ji ' 



(16) 

where k& is the Boltzmann constant and T is the tem- 
perature of the mechanical reservoir. The cavity field 
vacuum input noises have the following nonzero correla- 
tion functions 



k 

• k r (*K n (*')) 



in,t/ 



[N(u c ) + l]6{t - t'), (17a) 
N(u c )5(t-t'), (17b) 



with k = 1,2 assigning the two noise ports [see Fig. 1], 
and N(ui c ) = [exp{hui c / k#T} — 1] _1 is the mean intra- 
cavity photon number at thermal equilibrium with tem- 
perature T. At optical frequencies hw c /k-QT 3> 1 and 
therefore N{uj c ) w 0, so that only the correlation func- 
tion of Eq. (fTTa) is relevant. 

As mentioned above, our aim is to focus on the strong 
optomechanical regime (both radiatively and photothcr- 
mally) which establishes strong quantum correlations be- 
tween the steady state mechanical and optical fluctua- 
tions. It is shown in Refs. [19( and [20J that this is at- 
tained for very intense intracavity fields, i.e., for high 
finesse cavities and enough driving powers. In this case, 
when the system is stable, it reaches a steady state which 
is characterized by the cavity mode in a coherent state, 
and the mechanical resonator at a new equilibrium po- 
sition which is determined by the stationary intracav- 
ity photon number. Thus, it is enough to focus onto 



the dynamics of the quantum fluctuations around their 
classical steady state. This situation simplifies the exact 
QLEs in (fTS"]) to a set of linearized equations for quan- 
tum fluctuations of the system. Therefore, one replaces 
a generic operator O with O = O s + SO in Eq. (|15p , and 
gets the steady state values as a s = £/ {y c + iA) and 
q s = agGo(l + 272X/3£)/w m . Here the effective cavity 
detuning is introduced as A = Ao — Goq s , meanwhile we 
have chosen the phase of £ such that a s is real and pos- 
itive without loss of generality. Thus, dynamics of the 
quantum fluctuations are given by the following equa- 
tions 



Sq 
Sp 



u m Sp, 



(18a) 



G 



hit- s){2^ 2 Sx- 

7th 

— Lu m Sq — j m Sp + GSx + £, (18b) 
Si = -j c 5x + ASy + y^x™ + y/Zy^xf, (18c) 
Sy = - lc Sy - ASx + GSq + y/WiVi + y^vll^d) 

where the effective optomechanical coupling is defined as 
G = y/2 a s Go- In fact, G is a measure for optome- 
chanical coupling strength which is dominantly deter- 
mined by number of intracavity photons. Furthermore, 
in ([T5]) we have introduced the cavity field quadratures 
Sx = {5a + Sa^)/\/2 and Sy = (Sa — 5ai)/i\/2 and corre- 
sponding noise operators x™, y™ with k — 1,2 . 



B. Effective damping and optical spring effect 

Since the dynamics of the system is linearized and also 
the system noises are Gaussian, we can focus on a set 
of states; that is Gaussian states which are very inter- 
esting both experimentally and theoretically. In partic- 
ular, the bipartite steady state of the system formed 
by the mechanical oscillator and the cavity mode is 
a zero mean Gaussian state Hence, this state 

is fully characterized by its 4 x 4 correlation matrix 
(CM). Defining vector of the system operators as uit) = 
[Sq{t) 7 Sp(t),Sx{t),Sy(t)] T , the CM is given by 

Vij = -(uj(oo)itj(oo) +Uj(oo)uj(oo)}. (19) 

One may determine the CM of the system by solving 
Fourier transform of Eqs. (|18p . which are the following 
equations 

— iuoSq — w m <5p, (20a) 
—iujSp = —uj m Sq — 7 m <5p + GSx 

+V2a s \h(u) (272<5£ - y/VftS?) + £ (20b) 
-iwSx = -7 c <5x + ASy + \/2^x^ + y/2^x%, (20c) 
-iuSy = - lc 5y - ASx + GSq + yWiVi + vWW d ) 
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Subsequently, the CM elements can be obtained as 
da; da/ 



x 10 



)tl~ 



-in 



(ui(ui)Uj(uj') + Uj(u}')Ui(u)), 

(21) 

In particular, the first element of the CM which cor- 
responds to covariance of the mechanical position, (Sq 2 ), 
reads 



(Sq 2 



did 

2^ 



(22) 



where S qq (uj) is position spectrum of the mechanical res- 
onator, given by 

S qq {ui) = |X e ff(w)| 2 [5th(w) +S rp M +5pt(w) + <S cc (o;)], 

(23) 

where X c g(w) is the effective mechanical susceptibility. 
The spectrum of thermal Brownian force is (47l . l48l | 



<5 t h(w) 



j m ui fro; 
coth[ 



2fc B T J ' 



(24) 



while the radiation pressure and the photothermal spec- 
tra are given respectively by the following 



s rD M 



7cG 2 ( 7c 2 



- OJ 



[7 c 2 + (c-A)2][ 7 2 + (c + A)2] 
7 2 G 2 (x/3i) 2 



(25a) 
(25b) 



Finally, cross correlation of the radiation pressure and 
photothermal is introduced by the following spectrum 



S cc (uj) =2 l2 G 2 X f3L 

„ 7 c (7 2 + A 2 + ^ 2 )-W( 7c 2 



(26) 



■ U3 



(1 + T t V)[ 7 2 + ( W - A)2][ 7 2 + ( W + A)^] ■ 

In fact, the photothermal force in an optomechanical sys- 
tem -in addition to its exclusive effect- interacts with 
the radiation pressure. Such a cross correlation effect 
means that these two light forces do not act on the me- 
chanical resonator, independently. This is not surprising, 
since from (|18p one can find out that both radiative and 
photothermal forces are coupled to the mechanical res- 
onator through in phase quadrature of the cavity field 
Sx . However, comparing Eqs. (125b ) and (I26[) , the im- 
portance of the joint effect may be clarified. By investi- 
gation, one finds out that around the natural frequency 
of the mechanical resonator w m the crucial factor is fixL 
whose quantity for typical photothermal setups is much 
less than unity. Thus, at least around mechanical fre- 
quency the cross correlation effect has more contribu- 
tion than the pure photothermal spectrum. Nonetheless, 
around the mechanical frequency, the cross correlation 
spectrum is very sensitive to 71, which limits the value 
of the cross correlation spectrum. To visually show this, 
we compare the three spectra in Fig. 2, from which it is 
obvious that especially around the mechanical frequency, 




co/co 
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FIG. 2. (Color online) Mechanical Spectrum composition; 
radiation pressure (blue circle), photothermal force (red 
square), cross correlation effect (green plus). The system is 
composed of a mechanical oscillator with natural frequency 
cj m /27r = 10 MHz, quality factor Q m = 10 J , and mass m = 5 
ng. Temperature of the reservoir interacting with the me- 
chanical resonator is T = 0.4 K. Deformation constant of the 
cantilever is \ — 1 x 10 s/m, while its quantum efficiency 
is = 1 and the thermal diffusion time rth^m = 1. The 
cavity length is L = 1 mm and it is pumped by a laser with 
wavelength Al = 810 nm and power of P = 1 mW. The laser 
detuning from the corresponding cavity mode is A/o; m = 1. 
Finally, the optical damping rates are 71 /w m = 0.5 and 
72/71 = 1 • 



the cross correlation spectrum is even larger than the 
photothermal spectrum. One can also see from the fig- 
ure that destroying the total noise spectrum at high fre- 
quencies is the other consequence of the cross correlation 
effect. In Fig. 2 we have considered a set of parameters 
which are generally available in laboratories (see caption 
of Fig. 2). 

The effective mechanical susceptibility introduced in 
Eq. (gU) is defined as 



Xcff(w) = — 
0J m 



Sl 2 ff - - iuT e f{ 



(27) 



at which the effective resonance frequency f2 e ff and ef- 
fective mechanical damping rate r o g, are specified. This 
equation states that the mechanical damping rate is jus- 
tified by the cavity field forces. Splitting the effect of 
radiation pressure and bolometric forces, the effective 



damping rate reads as T 



bS 



In 



r rp + r pt , where 



2 7c AG 2 w n 



Fl ' p " [ 7c 2 + ( Wm - A)3] [ 7c 2 + (w m + A)']' (28a) 
2 l2 (3 x LAG 2 uj m [2 lc + r th (A 2 + 7 2 - <)} 
pt " (1 + r t 2 h O[7 2 + ( Wm - A)2][ 7 2 + ( Wm + A?* 



|b) 



On the other hand, the modified mechanical frequency, 
which is induced by light (optical spring effect), is given 
by 



7 



^eff = 



AG 2 Wm ( 7c 2 + A 2 - w 2 ) 



2 72 /3 X LAG 2 o; m [7 c 2 + A 2 - lj 2 (1 + 2 7c r th )] 



[ 7c 2 + K 1 ~A) 2 ][ 7c 2 + K 1 + A) 2 ] (! + • 



th^m 



)[ 7c 2 + K 1 -A) 2 ][ 7c 2 + K 1 + A) 2 



(29) 



This can lead to significant change in the mechanical fre- 
quency only close to bistability threshold, but as far as 
the system works at weak coupling regime modification of 
the mechanical frequency is negligible. It is straightfor- 
ward to verify that Eqs. ([28]) and (|29f reproduce previous 
results [39|,[4(|. Moreover, by setting 7 2 equal to zero, the 
effective frequency and effective damping rate of a pure 
radiation pressure optomechanical system are obtained 
which are already reported in Ref. [37| . 



IV. PHOTOTHERMAL COOLING 

In this section, we study cooling the motion of the 
mechanical resonator toward its ground state by means 
of the correlation matrix introduced in Sec. III. After- 
wards, to have a physical intuitive picture and to com- 
pare our results to the other works, we consider a weakly 
coupled optomechanical system to approximate the exact 
results. However, before every thing, the stability condi- 
tions must be carefully taken into account. Because our 
steady state analysis is valid only for a stable system. 

The system is stable when all poles of the effective sus- 
ceptibility ([27|) lie in the lower complex half plane. The 
parameter region in which the system is satisfactorily sta- 
ble can be obtained from the Routh-Hurwitz criteria [4!| . 
The inequalities determining the stability conditions are 
too cumbersome to be reported here. However, one can 
appreciate this from the characteristic polynomial of the 
effective susceptibility (see Appendix B). In the case of 
a red detuned cavity (A > 0), which is the interesting 
regime both for cooling and entanglement, the only non- 
trivial stability condition is 



tion variance by 



A 2 ) 



AG 2 (l + 2 72 f3 x L) > 0. 



(30) 



In the rest of the paper satisfaction of this stability con- 
dition is considered. 

The effective phonon occupation number of the me- 
chanical resonator, n e ff, can be derived from the total 
mechanical energy 



= —frbJ, 

2 



.[(dp 2 ) + (Sq 2 )] = hw m {n cif + ~). (31) 



The mechanical resonator approaches to its ground state 
if n e ff ~ or equivalently E m ~ tkJ m /2. Therefore, to an- 
alyze the cooling process, one needs to calculate the me- 
chanical position and momentum variances, which are re- 
spectively first and second diagonal elements of the CM. 
We have already discussed about (Sq 2 ), while variance of 
the mechanical momentum is simply related to the posi- 



(SP 2 ) = I ~ S «M- 



(32) 



The integrals in (|22p and (|3"2"|) can be performed for a 
stable system. However, a sensible approximation on the 
Brownian noise equation, significantly simplifies the final 
result. In fact, the Brownian noise £(i) is in general a 
non Markovian Gaussian noise [see Eq. (|24p] . Neverthe- 
less, to observe the quantum mechanical effects of the 
mechanical resonator, it is necessary to prepare a high 
quality mechanical resonator Q m 3> 1. At this limit 
becomes Markovian and we have 



coth[— — ] 



7n 



2fc B r 

1 fajJm 



7 m (2n + l) 



(33) 



where n = [exp(/iu; m /fcBT)~l]~ 1 is mean thermal phonon 
number of the mechanical resonator connected to a reser- 
voir at equilibrium temperature T. 

Another important point must be considered here is 
that the ground state cooling is not approached only by 
attaining n e ff < 1 . In fact, in order to get the quan- 
tum ground state cooling, it is also necessary for both 
the mechanical displacement and momentum variances 
to satisfy the minimum uncertainty principle. Strictly 
speaking, both variances have to tend to 1/2, therefore, 
energy equipartition has to be satisfied in the optimal 
regime close to the ground state. It is essential to be 
pointed out that, in general, energy equipartition is not 
necessarily true as it can be seen from (f2"2|) and (f3"2"]l that 
the two variances are not equal (Sq 2 ) ^ (Sp 2 ). 



A. Weak coupling limit 

Performing the integrals for (Sq 2 ) and (Sp 2 ), and then 
replacing in (|3"Tj) . gives an exact result for the effec- 
tive mechanical phonon number, which is too involved 
to be reported here. However, at weak coupling regime 
7c 3> 7m, G, one gets an approximate expression for the 
effective phonon number. Additionally, the cooling pro- 
cess can be executed at low temperatures, at which num- 
ber of the thermal phonons is very small w m 3> n"f m - 
Exerting the above limits; weak optomechanical regime 
and working at cryogenic temperatures, we arrive at the 
following approximate equation for the effective phonon 
number 



"■off 



n 7m + A+ + A+ t 



7m + r rp 



pt 



(34) 
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where we have introduced the photon scattering rates 
A+, Ap t , and A+ corresponding to the radiation pres- 
sure, photothermal force, and their cross correlation ef- 
fect, receptively. These are the rates at which the laser 
photons are scattered by the mechanical oscillator simul- 
taneously with the absorption of the oscillator vibrational 
phonons, given by 



a: 



g 2 1c 



2[7 C 2 



p* ~ 2(1 



r t 2 h W m 
2 



A nr . — 



7 2 /3x^G' 2 [7 c - T th ^ m (A + u>„ 
(l + ^<)[7c 2 + K + A) 2 



(35a) 
(35b) 
(35c) 



Interestingly, it is obvious from Eqs. (|35b ) that for a red 
detuned cavity, depending on the working regime, A^ c 
can get positive or negative values. Actually, for an op- 
tomechanical system in bad cavity regime ( 7c 3> w m ) 
gets positive values, while for a good cavity system 
( 7c -c u> m ) it becomes negative. Such a change in the 
sign of ^4+ at the good cavity regime can play an im- 
portant role in the final effective phonon number. It is 
notable that in this regime the radiation pressure cools 
down the mechanical resonator very effectively [35l [3(| . 

One can also obtain Eq. (|34p by employing the asym- 
metric noise spectrum method exploited for the radiation 
pressure cooling in Ref. 35]. According to the method, 
the minimum attainable phonon number is given by the 
following equation 



(+w m ) 



(36) 



where S op t(u)) is the asymmetric noise spectrum of the 
optical force. In Ref. [4JJ the noise spectrum of the light 
forces, including both radiation pressure and photother- 
mal pressure, was applied to Eq. (f36| to get the effective 
mechanical phonon occupation. Considering the differ- 
ence in definition of the elasticity constant (one needs 
to replace 2/3/c by 0x), coincidence of our result with 
the result of Ref. [4l| can be easily verified. Moreover, 
it can be shown that suppressing A+ c in Eq. (jM]) gives 
the result obtained in Ref. 40]. In fact, in Ref. [4(| it 
was supposed that the radiation pressure and photother- 
mal force operate independently, so the cross correlation 
effect in photothermal cooling is omitted. Our solution 
is also in coincidence with the quantum theory of radi- 
ation pressure cooling [35T - I371 ] . Actually, if no photon is 
absorbed by the micromechanical mirror, i.e. 7 2 = 0, 
Eq. (|3"4")l reproduces the result obtained by Genes et al. 



B. Cooling at strong coupling regime 

Now let us go beyond the previous works and study 
the photothermal cooling at strong coupling regime. We 
use the exact effective phonon number obtained from 
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FIG. 3. (Color online) Variation of the effective phonon 
number versus normalized cavity detuning for a bad cavity 
7i/cj m = 10 (a) and a good cavity 7i/o; m = 0.1 (b): Bare radi- 
ation pressure system (blue circle) and for a system including 
both radiation pressure and photothermal effects; 72/71 =0.1 
(green plus), 72/71 = 0.5 (red square) and 72/71 = 1.0 (cyan 
triangle). The contour plot of inset (c) presents the effective 
phonon number versus normalized thermal diffusion time and 
normalized decay rate from the micromechanical mirror side 
(photon absorption rate) for a good cavity (71 /w m = 0.1) at 
fixed cavity detuning A = w m . The input power is P = 15 
mW and the rest of the system parameters are same as the 
parameters used in Fig. 2. 



Eq. (|3"Tj) to make a numerical investigation on the cooling 
process. In the following analysis, we implicitly inspect 
satisfaction of the energy equipartition which is, as men- 
tioned above, an important condition for ground state 
cooling. In fact, in the parameter region utilizing in this 
paper -both for bad and good cavities- it is always true 
that (Sq 2 ) ~ (Sp 2 ). 

In Fig. 3(a) and (b), effective phonon number of the 
micromechanical resonator is plotted versus cavity de- 
tuning for a bad cavity and a good cavity, respectively. 
The parameters we are considering here correspond to 
a practically attainable setup, for example one can op- 
erate the mechanical resonator in Ref. |5G ] as the mov- 
able mirror of an optomechanical system similar to the 
setup in Ref. [5l|. The figure shows that for a bad cav- 
ity the minimum achievable phonon number of a pho- 
tothermal system is much smaller than a pure radiative 
system. The minimum number of phonons in this case 
is about five phonons which is exactly the same as the 
result of Ref. [40]. On the other hand, surprisingly the 
photothermal force has a positive effect on the ground 
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state cooling even in a good cavity setup. This can be 
found out from Fig. 3(b), where the minimum phonon 
number accomplished by a pure radiation pressure sys- 
tem is about 0.086, while adjustment of photothermal 
effects (by adjusting the rate of photon absorption at 
the micromechanical mirror) reduces it to about 0.075. 
Therefore, better cooling results can be achieved by engi- 
neering the photon absorption rate of the micromechan- 
ical mirror. For example, one can control the absorption 
rate 72 by spreading gold nanoparticles on the surface of 
the cantilever. Adjusting the nanoparticle size, justifies 
its optical absorption coefficient, which in turn modifies 
absorption rate of the mechanical resonator. 

However, the optomechanical cooling can be optimized 
more by another characteristic parameter of the pho- 
tothermal force, that is the thermal diffusion time r t h, 
which also functions as the delay time in the feedback 
loop. In Fig. 3(c) we have depicted variation of the ef- 
fective phonon number with the absorption rate, 72, and 
thermal diffusion time, Tth, for a good cavity system. It is 
evident from the figure that there is an optimum value for 
both absorption rate and thermal diffusion time, which 
gives the minimum effective phonon number. Employing 
these optimal values in a photothermal system, the effec- 
tive phonon number can even be reduced to 0.072. This 
is equivalent to an effective temperature about 178 fiK 
for the mechanical mode. 

It is also interesting to compare the exact result, which 
enables us to consider the system at strong coupling 
regime, with the approximate calculation. For the op- 
tomechanical parameters employed in this paper, a bad 
cavity system is always working at weak coupling limit. 
Therefore, the approximate result reproduces just the 
exact solution. However, for a good cavity, the weak 
coupling limit is only obtained for small input powers. 
Thus, we focus on the good cavity regime to see how can 
strong optomechanical coupling affect the ground state 
cooling. For this purpose, in Fig. 4 we have plotted the 
effective mechanical phonon number versus input laser 
power which mostly determines value of the effective op- 
tomechanical coupling. In the figure, the exact and ap- 
proximate values of n G ff are included both for a pure ra- 
diative system and a photothermal scheme. The dashed 
vertical green line indicates the laser power which gives 
G = 7c, therefore, the right hand side of the line ab- 
solutely refers to the strong coupling region. It can be 
seen that for P < 2 mW, which pertains to the weak 
regime, the exact and approximate results begin to coin- 
cide, while for larger values of P the approximate lines 
go far from their corresponding exact lines. Neverthe- 
less, one can conclude from Fig. 4 that the ground sate 
cooling is enhanced at strong coupling regime, especially 
photothermal effects are more responsive at this regime. 
Additionally, the figure shows that there is a finite opti- 
mal value for G which gives the minimum n e ff. Hence, 
an arbitrary large G does not necessarily improve the 
ground state cooling. 

Finally, let us check validity of our theoretical formu- 
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FIG. 4. (Color online) Variation of the effective phonon num- 
ber versus laser input power for a good cavity ji/ui m = 0.1. 
For a pure radiative system; exact solution (blue circle) and 
approximate solution (cyan filled circles). For a photother- 
mal scheme with 72/71 = 1; exact solution (red square) and 
approximate solution (orange filled square). The cavity de- 
tuning is A = u)m and the rest of the system parameters are 
same as the parameters used in Fig. 2. 

lation with already existing experimental evidence, and 
then investigate availability of better results by the same 
experimental setup. According to our results, the effec- 
tive temperature of the mechanical resonator is about 30 
K for 7i/w m = 800 and A/w m = 650. These are nearly 
the same experimental parameters utilized in Ref. [3Q| to 
achieve a cantilever cooled down to 30 K, therefore, our 
analysis gives a good insight for operation of the pho- 
tothermal effect. On the other hand, for the parameters 
of the setup in Ref. [2!| (very similar to that of Ref. [30| ) , 
we have computed effective temperature of the mechani- 
cal resonator for a wide range of cavity detuning A and 
the cavity input side decay rate 71. Our numerical re- 
sults indicate that, even though the authors in Ref. [2^] 
— starting from room temperature- achieved a minimum 
temperature about 70 K for the mechanical resonator, 
one is able to cool it down to effective temperatures as 
less as 0.2 K. In fact, employing the same device, it is 
possible to achieve such a low temperature by improving 
the cavity finesse such that 7i/u; m 60, and adjusting 
the laser frequency detuning. 



V. PHOTOTHERMAL ENTANGLEMENT 

Eventually, in this section we shall discuss about the 
most important result of the paper; the optomechanical 
entanglement in a system coupled through both radia- 
tion pressure and photothermal force. Entanglement be- 
tween various constituents of a radiation pressure based 
optomechanical system is thoroughly studied, including 
optomechanical entanglement of a mechanical resonator 
and optical cavity modes and also all optical entangle- 
ment of the output modes (see e.g. Ref. 20]). However, 
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effect of photothermal pressure on the optomechanical 
entanglement which can be an interesting concept is not 
studied, so far. The optomechanical entanglement is pro- 
duced only in a strongly coupled system, therefore, only 
a considerable contribution of the photothermal force can 
enhance it. In fact, although the photothermal force nat- 
urally introduces decoherence, yet a strong photothermal 
action even can improve the optomechanical entangle- 
ment. To investigate the situation, we first provide an 
exact numerical description, and then an approximate 
analytical analysis on the entanglement is presented. 

In the following we use logarithmic negativity as a con- 
venient measure for continuous variable entanglement. 
Log negativity is defined as [HI, [53[ 



Em — max[0, — In 2?? ], 



(37) 



where rf~ is the minimum symplectic eigenvalue of partial 
transposed CM corresponding to the bipartite system of 
interest, which is given by 

rf = 2- 1 / 2 {£(V0 - ^Z(V) 2 -4detV} 1/2 , (38) 

with E(V) = dctV^A + dct V B - 2det V C - Here we have 
considered the CM in a 2 x 2 block form as 



V 



V A V c 



(39) 



Therefore, from (|37|) and (|38|) one concludes that a Gaus- 
sian state is entangled if and only if rf~ < 1/2, which is 
equivalent to Simons necessary and sufficient entangle- 
ment nonp ositive partial transpose criterion for Gaussian 
states |54|. 



A. Exact numerical results 

Now let us study behaviour of entanglement between 
the mechanical resonator and the optical cavity mode, 
Em, from exactly calculated CM (IT91 . We shall focus on 
the two exclusive photothermal parameters, the charac- 
teristic time of the thermal diffusion r t h and the pho- 
ton absorption rate of the movable mirror 72, to ex- 
amine their effect on the optomechanical entanglement. 
We have made a careful analysis in a wide parameter 
range and we found a parameter region very close to per- 
formable photothermal experiments (Hoi . l5lj . 

In Fig. 5(a), variation of Em versus 72 and rth is shown 
by a contour plot. The lower boundary of the plot, cor- 
responding to 72/71 = 0, refers to a system which is 
coupled only by the radiation pressure. It is clear from 
the figure that for a radiative system, the log negativity 
is about 0.4, which -as expected- is independent of rth- 
At the first sight one can see that, interestingly, increas- 
ing 72, which brings in photothermal effects, can improve 
the entanglement. One can see from Fig. 5(a) that for 
small values of the thermal diffusion time (rthW m < 1), 
Em decreases rapidly with 72, therefore, it can be con- 
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FIG. 5. (Color online) The contour plot of logarithmic neg- 
ativity versus normalized thermal diffusion time and normal- 
ized decay rate from the micromechanical mirror side (a). 
Variation of log negativity versus normalized input mirror 
decay rate (b). Variation of log negativity versus tempera- 
ture of the mechanical reservoir (c): In insets (b) and (c) blue 
circle corresponds to a bare radiation pressure system, while 
the other lines belong to a system including both radiation 
pressure and photothermal effects; 72/71 = 0.2 (green plus), 
72/71 = 0.9 (red square). In the case of inset (a) and (c) the 
input side decay rate and the laser power are 7i/w m =0.1 
and P = 50 mW, respectively. In the inset (b), the cyan re- 
gion indicates the unstable region of a pure radiation pressure 
system and the input power is P = 20 mW. The cavity detun- 
ing is always fixed at A/tj m = 0.85. The rest of the system 
parameters are same as the parameters used in Fig. 2. 



eluded that a fast photothermal action -compared to the 
mechanical response time 1 / w m - severely destroys the en- 
tanglement. However, at larger rth values, Em reaches to 
its maximum for a photothermal system (72 7^ 0). This 
consequence means that in some cases, the photothermal 
force can lead to enhancement of the optomechanical en- 
tanglement. Especially, in the parameter regime selected 
in this paper, the log negativity increases significantly to 
0.5 for an absorption rate 72 equal to the input mirror 
decay rate 71 . 

Instead, in Fig. 5(b), we plot the log negativity re- 
spect to the input decay rate 71 at fixed pump power 
and cavity detuning. The teal lines indicate the param- 
eter region at which the pure radiation pressure system 
is unstable. Interestingly, from the figure it is obvious 
that the photothermal force can significantly modify the 
stability region. In fact, the optomechanical entangle- 
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ment of a pure radiative system cannot raise enough due 
to instability of the system, while in the presence of pho- 
tothermal effects the bistability threshold retreats, there- 
fore, En can raise to its maximum. However, the raising 
slope of En in the photothermal case is smaller than the 
radiation pressure case, since the photothermal force in- 
troduces decoherence to the system. According to the 
figure, further increase in 72 leads to destruction of the 
entanglement, as decoherence of the system grows with 
it. 



Now let us investigate the interesting concept of the 
entanglement robustness against temperature of the me- 
chanical reservoir. To show effect of the photothermal 
force on the robustness, in Fig. 5(c) the log negativity 
of radiation pressure and two different photothermal sys- 
tems is depicted versus temperature of the mechanical 
resonator environment. The figure shows that in addi- 
tion to the amount, introducing an absorbing microme- 
chanical mirror to the optomechanical system can even 
increase robustness of En- Thus, one can enhance the 
robustness of En by controlling absorption rate of the 
micromechanical mirror 72, which can be achieved by 
spreading metallic nanoparticles on the surface of the 
mirror. 



B. Approximate analytical results 

The exact expression of the logarithmic negativity 
stemming from the solution of (|21j) is very cumbersome, 
but it is nonetheless possible to explain some of the above 
results by means of an approximate treatment which sat- 
isfactorily describes the effect of photothermal force on 
En- We derive an approximate analytical expression for 
En in the parameter region very close to the bistabil- 
ity threshold. This regime is relevant for optomechanical 
entanglement because, as first pointed out in Ref. [2(| 
and later discussed in detail in Ref. (55|, En reaches 
its maximum value at the bistability threshold. There- 
fore, taking the threshold value for the optome chanical 
coupling G ~ v/w m (7 c 2 + A 2 )/[A(1 + 2 l2 f3 X L)], and ne- 
glecting thermal noise terms one gets the following ap- 
proximate expression for the minimum symplectic eigen- 
value 

_ r 7cK + 4A 2 (A 2 + 7c 2 + Q] A v + p X L>y 2 B v y h 

11 I 2[87cA 2 (A 2 +7 c 2 + 5u; 2 n )A, ; + ^72C, ) ] J ' 

(40) 

where 

A v = l + rth [A 2 T th + 7 C (2 + 7c 7th) + w 2 1 r th (l + 2 7c Tth)] , 
B v = b + b 2 0J^ n + biUJ^ + beu;^, 
C v = 2b + c 2 o4 + c 4 u4, 

in which we have 



60 =40A 2 7c (A 2 + 7c 2 )(l + 2 7c r th + A 2 r 2 h + 7c 2 r 2 h ), 

b 2 =8A 2 [6 7c 4 r t 3 h - A 2 r th + 7c 5 r t 4 h + 7c 3 r 2 h (3 + 2A 2 r 2 h ) + 7 2 Tth(7 + 6A 2 r 2 h ) + 7c (5 + 7A 2 r 2 h + A 4 r t 4 h )] , 
h =2 [ 7c 3 T 2 h (3 + 4A 2 r 2 h ) + 8 7 2 (r th + 3A 2 r t 3 h ) + 7c (5 + 7A 2 r 2 h + 4A 4 r t 4 h ) - 6A 2 r th ] , 
b 6 =2r th [6 7c 2 r 2 h + 7 c 3 T t 3 h + 7cT t h(A 2 T 2 h - 1) - 2] , 

c 2 =16A 2 [I0 7c 4 r 3 h + 7 c 5 r 4 h + 7c 3 T 2 h (27 + 2A 2 r 2 h ) + A 2 r th (3 + 4A 2 r 2 h ) + 7 c 2 r t h(43 + 14A 2 r 2 h ) 

+ 7c(21 + 3lA 2 r 2 h + A 4 r 4 h )], 
c 4 =16A 2 r th [2A 2 r 2 h + 247 c 2 r 2 h + 5 7c 3 r t 3 h + 7crt h (5A 2 r 2 h - 3) - 7] . 



If the photothermal effect is negligible, 72 
gets 



+ 4A 2 ( 7c 2 +A 2 



0, one 



(41) 



16A 2 (A 2 +7c 2 + 5<)] 



which is in coincidence with the results of Refs. |55l l56j. 
Instead, in the presence of the photothermal force the in- 
equality 72 < 7c imposes that the terms proportional to 
72 in Eq. (|40|) cannot make a predominant contribution 
either in the numerator and denominator. Hence, accord- 
ing to the parameter set employed in this paper, the pho- 



tothermal force can slightly modify the optomechanical 
entanglement produced by the radiation pressure. Such 
a modification in some cases can lead to enhancement 
of the entanglement. In fact, bo is the prevailing factor 
both in B n and C v , while in the latter it is multiplied by 
two. Therefore, in general it is true that B n < C n , which 
means that the value of r)~ becomes less (so En becomes 
greater) than its value for a pure radiation pressure sys- 
tem. However, depending on the system parameters it 
is also possible to have B v > C v , and therefore rj~ may 
easily tend to 1/2 (so En tends to zero), even for not 
too large values of the absorption rate 72. Nevertheless, 
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the other effect of the photothermal force, as discussed 
above, is adjustment of the bistability threshold which 
allows for raising of the En, meanwhile a pure radiative 
system is unstable. Thus, the entanglement in an pho- 
tothermal system can be higher than a mere radiation 
pressure scheme even if B v > C v . 



VI. CONCLUSION 

We have studied the entanglement of a micromechan- 
ical mirror with a cavity mode and cooling of the res- 
onator in the presence of photothermal effects. We have 
analyzed the dynamics by adopting a standard quantum 
Langevin treatment. Dynamics of the quantum fluctua- 
tions which is expressed by the linearized QLEs around 
the classical stationary state of the system have been an- 
alyzed. We have investigated dependence of the log nega- 
tivity and the mechanical occupancy of the optomechani- 
cal stationary state upon various system parameters. We 
have also derived approximate, but compact, analytical 
expressions for effective mechanical phonon number in 
a photothermal system at weak coupling regime and for 
the minimum symplectic eigenvalue of the bipartite sys- 
tem (composed of mechanical resonator and cavity mode) 
close to the bistability threshold. At the strong coupling 
regime we have shown numerically that the cooling is im- 
proved by the photothermal coupling. Finally, it is real- 
ized from a numerical study and verified analytically that 
even though bolometric force is naturally a dissipative ef- 
fect, it can improve the optomechanical entanglement. 
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Appendix A 

The thermoelastic wave equation is given by 
Ey _ . , Ey 



gu 



"2(l + cr)(l - 2a) 
E Y a t h 



V(V.u) 



2(1 + 0") 



V 2 u 



3(1 - 2a) 



VT, 



(Al) 



where u is the displacement vector, Ey is Young mod- 
ulus and a is Poisson ratio of the material construct- 
ing the mirror [42 ] . Splitting the displacement vector 
into transversal and longitudinal components, that is 
u = u; + u t where V.Uf = and V x u/ = 0, the wave 
equation for every component reads 



3g(l - 2a) 



u t - c 2 V 2 u t = 0, 



where 



ct 



Ct 



, E Y (l-a) 
[ g(l + a)(l-2a)' 

r E Y il/2 

l 2g(l + aY ' 
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(A2a) 
(A2b) 

(A3a) 
(A3b) 



are speed of longitudinal and transverse thermoelastic 
wave propagation, respectively. 



Appendix B 



Characteristic polynomial of the effective susceptibility 
Eq.fl27|) is P(tp) — ip 5 + pxip 4 + p 2 ip 3 + P3<p 2 + Pif + P5, 
where 



Pi = 7m + 27 c + 



1 

Tth 



(Bla) 



P2 = A 2 + < + 7c 2 + 2 7c7m + 2iH±r> , (Bib) 



Tth 



p 3 = — {A 2 (l + 7m r th ) 



+7c[7c + 7m (2 + 7cT th ) + 2T th o;i]}, (Blc) 



P4 = — {(7m + TthW^XA 2 + 7c ) 

Tth 

+w m (2 7c u; m - G 2 Ar th )}, 



(Bid) 



p 5 = — {<(A 2 + 7c 2 ) - G 2 A^ m (l + 2 72X /3L)}Blc) 

Tth 
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